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ABSTRACT 

We continue our study of weakly ionized protostellar accretion discs that are threaded 
by a large-scale magnetic field and power a centrifugally driven wind. It has been 
argued that there is already evidence in several protostellar systems that such a wind 
transports a significant fraction of the angular momentum from at least some part of 
the disc. We model this situation by considering a radially localized disc model in which 
the matter is everywhere well coupled to the field and the wind is the main repository 
of excess angular momentum. We consider stationary configurations in which magnetic 
diffusivity counters the shearing and advection of the magnetic field lines. In Wardle 
& Konigl we analysed the disc structure in the hydrostatic approximation (vertical 
motions neglected inside the disc) and presented exact disc/wind solutions for the am- 
bipolar diffusivity regime. In Konigl, Salmeron & Wardle (Paper I) we generalized the 
hydrostatic analysis to the Hall and Ohm diffusivity domains and used it to identify 
the disc parameter sub-regimes in which viable solutions with distinct physical prop- 
erties can be expected to occur. In this paper we test the results of Paper I by deriving 
full numerical solutions (integrated through the sonic critical surface) of the disc equa- 
tions in the Hall domain. We confirm all the predictions of the hydrostatic analysis 
and demonstrate its usefulness for clarifying the behaviour of the derived solutions. 
We further show that the outflow solutions can be extended to larger scales (so that, 
in particular, they also cross the Alfven critical surface) by matching the localized 
disc solutions to global 'cold' wind solutions of the type introduced by Blandford & 
Payne. To facilitate this matching, we construct a library of wind solutions for a wide 
range of wind model parameters; this library is made available to the community. 

The results presented in Wardle & Konigl, Paper I and this work combine to form 
a comprehensive framework for the study of wind-driving accretion discs in proto- 
stellar and other astrophysical environments. This theoretical tool could be useful for 
interpreting observations and for guiding numerical simulations of such systems. 

Key words: accretion, accretion discs - ISM: jets and outflows - MHD - stars: 
formation. 



1 INTRODUCTION 

A common feature of protostellar accretion discs is their 
association with collimatcd, energetic outflows (e.g. Cabrit 
2007). These outflows are widely believed to represent cen- 
trifugally driven winds that are launched along large-scale, 
ordered magnetic fields (e.g. Konigl & Pudritz 2000; Pu- 
dritz et al. 2007). Such winds could in principle be efficient 
trasnporters of disc angular momentum (e.g. Blandford & 
Payne 1982, hereafter BP82), and it has in fact been ar- 
gued that there is already observational evidence in several 
protostellar systems that an outflow of this type carries the 



bulk of the excess angular momentum from at least some 
part of the associated accretion disc (e.g. Ray et al. 2007). 
Protostellar discs are typically weakly ionized, and a certain 
minimum degree of ionization is required to attain the level 
of coupling between the matter and the field that enables the 
vertical magnetic angular-momentum transport mechanism 
to operate. [A similar requirement must be satisfied also to 
enable radial angular-momentum transport by a small-scale, 
turbulent magnetic field; such turbulence could be induced, 
for example, by the magnetorotational instability (MRI; e.g. 
Balbus & Hawley 1998).] The inherent magnetic diffusivity 
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tends to counter the efTects of shearing and advection by the 
differentially rotating accretion flow and therefore makes it 
possible for the discs to attain a steady state, at least on the 
dynamical (rotation) time. 

The behaviour of a weakly ionized gas can be charac- 
terized according to the nature of the dominant diffusivity 
mechanism for the given density and ionization state: am- 
bipolar, Hall, or Ohm (e.g. Konigl & Salmeron 2011). The 
ambipolar regime, in which the magnetic field lines are effec- 
tively frozen into the ions and drift with them relative to the 
dominant neutral component, can be expected to dominate 
over the entire disc cross section in the outermost regions of 
the disc (at radii r > 10 au) and near the disc surfaces at 
smaller radii. In the Hall regime the magnetic field is frozen 
into the electrons and drifts with them relative to the ions 
and neutrals; this regime could dominate over most of the 
disc cross sections on scales r ~ 1 — 10 au. The Ohm regime, 
in which the field lines completely decouple from the charge 
carriers, could potentially dominate near the disc mid-plane 
on scales ~ 0.1 — 1 au. (At smaller radii the mid-plane re- 
gion likely becomes coUisionally ionized.) However, as wind- 
driving discs typically have comparatively lower column den- 
sities and mid-plane densities because of their high angular 
momentum transport efficiency (which results in relatively 
high inflow speeds), the low- ionization Ohm regime might 
be curtailed (or even entirely eliminated) in such systems. 
Even if this were to happen, the Hall-dominated zone would 
be unlikely to extend to much smaller radii than estimated 
above on account of the fact that, at sufficiently large den- 
sities, the dominant charge carriers become positively and 
negatively charged grains of equal mass (e.g. Nishi, Nakano 
& Umebayashi 1991), so that no Hall current can fiow. 

The structure of radially localized wind-driving proto- 
stellar discs in the ambipolar diffusion-dominated regime 
was investigated by Wardle & Konigl (1993, hereafter 
WK93). They derived exact solutions that were matched to 
global, radially self-similar disc-wind solutions of the type 
introduced by BP82. They analysed the disc solutions by 
using the hydrostatic approximation, in which the vertical 
velocity component is neglected inside the disc, and obtained 
useful algebraic relations that led to a set of parameter con- 
straints on physically viable configurations (see also Konigl 
1997). This analysis was generalized by Konigl, Salmeron 
& Wardle (2010, hereafter Paper I) to the Hall and Ohm 
diffusivity regimes. In particular, they found that all the vi- 
able solutions correspond to four sub-regimes in the Hall 
domain and three in the Ohm domain, with the solutions 
in each of the identified parameter sectors having distinct 
physical properties. The main goal of this paper is to derive 
exact wind-driving disc solutions that can test these predic- 
tions and, more generally, the applicability of the hydrostatic 
approximation to the description of the salient features of 
such systems. We concentrate on the Hall diffusivity domain, 
which, according to the discussion above, should be most 
relevant (together with the ambipolar regime already dis- 
cussed in WK93) to the study of the weakly ionized regions 
of wind-driving protostellar discs. In Paper I we formulated 
the problem in terms of a conductivity tensor, specified by 
the values of the Pedersen, Hall and Ohm components, and 
discussed its correspondence to a multifiuid formulation for 
the case where only two charged species (one positive and 
one negative) are present. Our parameter constraints were, 



in fact, derived in the context of the latter approach. In 
this paper we employ the conductivity-tensor formulation 
for the characterization of the solutions that we derive, but 
we again employ the multifluid formulation in the analysis. 

The paper is organized as follows. Section 2 summarizes 
the system of equations that underlies the radially localized 
disc model and describes the model parameters as well as 
the boundary conditions on the equations and the method 
of their numerical integration. Section 3 provides the corre- 
sponding description of the global, self-similar wind model 
and outlines the procedure used to match the disc and wind 
solutions. Section 4 presents representative solutions for Hall 
diffusivity-dominated discs in the different parameter-space 
sub-regimes identified by the hydrostatic analysis in Paper I 
and compares them with the predictions (summarized in Ap- 
pendix A) of the analytical treatment. The main findings of 
the paper are recapitulated in Section 5. 



2 RADIALLY LOCALIZED DISC MODELS 

2.1 Dimensionless system of equations in z 

In Section 1.3 (where the prefix T hereafter denotes Paper I) 
we derived the algebraic relations and ordinary differential 
equations (ODEs) in the vertical cylindrical coordinate z 
that describe the vertical structure of the disc at any given 
value of the radial cylindrical radius r. For clarity and ease of 
reference, we reproduce them below in dimensionless form. 
They are the equations of motion 
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the relations linking the electric field in the inertial coordi- 
nate system and in the frame comoving with the neutrals 



and Ohm's Law 



Cr = WEr + Wcf, — WzO^, , 
= -CB + Wz&r - Wr ; 

CTH / / , 



(7) 
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jr = 2/(5-0 - ^p)br + — {e'zb^ - e'^) -|- crpe^ , (9) 



j<l> = Vi^o - a-p)60 + ^ (fir - e.'zbr) + ffve'^ , (10) 
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, _ —{e'rbr + e'^h^){ao — ffp) 
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where we have taken jz to be ~ 0. The following notation 
has been used in the above expressions: 
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where p is the gas density, v is the fluid velocity, J is the cur- 
rent density, B is the magnetic field, <r is the conductivity 
tensor, incorporating the Pedersen, Hall and Ohm compo- 
nents (cTp, cth and ao, respectively), and 



E' = E + 



V X B 



(15) 



is the electric field in the frame of the neutrals, in terms of 
E, the electric field in the inertial ('laboratory') frame. Also, 
y = E' -B/B^, vk is the Keplerian speed, Cg is the isothermal 
speed of sound, /it = (f's/cK)?' is the tidal scaleheight and 
the subscript '0' denotes the disc mid-plane. 

Under the thin-disc approximation, the vertical com- 
ponent of the magnetic field is constant with height; there 
is therefore no need to solve the vortical component of 
Ampere's Law. It can also be shown that the azimuthal com- 
ponent of E (or, equivalently, the variable WEr) is constant 
with height in our model (see Section 1.3.7). Therefore the 
variable vbt ~ cE^/Bz, which represents the radial drift 
velocity of the poloidal flux surfaces (see Section 1.3.6) is 
also constant with height inside the disc. 

In addition to the above equations, the following are 
also included so as to enable the position of the sonic point 
(subscript 's') and the (normalized) upward mass flux Wzo 
(= Ps under the assumptions that pWz = const and that the 
disc is vertically isothermal) to be derived self-consistently 
as part of the solution (see Section 2.4): 

0, (16) 



dz 
dz 
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2.2 Parameters 



The solutions are specified by the parameters enumerated 
below. Only a brief summary is provided here; the reader is 
referred to Section 1.3.13 for further details. 

(i) tto = vazo/cs, the mid-plane ratio of the Alfven speed 
(based on the uniform vertical field component Bz) to the 
isothermal sound speed. It is a measure of the magnetic field 
strength. 

(ii) Cs/vk ~ hT/r, the ratio of the tidal scaleheight to 
the disc radius - a measure of the geometric thinness of the 
disc. Although this parameter does not appear explicitly in 
the normalized equations, it is required for mathcing the 



disc and wind solutions (see Sections 3.2 and 3.5) and is 
used to place an upper limit on the midplane radial speed 
in physically viable solutions (see Appendix ??). 

(iii) The mid-plane ratios of the conductivity tensor com- 
ponents: ((Tp/crj_)() (or ((th/o"x)o) and {u±/ao)o, where 
a± = (an + o"p)^^^ is the total conductivity perpendicular 
to the magnetic field. They characterize the conductivity 
regime of the fluid and are taken here to be independent 
of 2 to facilitate the comparison with the analytic results 
derived in Paper I. 

(iv) The mid-plane value of the Elsasser number Ao = 
VAo/{^Kri±o), where r]±o = c^/4-!ra±o is the 'perpendicular' 
magnetic diffusivity and f^K ~ vk/t is the Keplerian angu- 
lar velocity. This parameter measures the degree of coupling 
between the neutrals and the magnetic field, with values ^ 1 
and -C 1 corresponding to strong and weak coupling, respec- 
tively. In this paper we assume (as was done in paper I) that 
the disc is everywhere magnetically well coupled, so that Ao 
is never <C 1. As an alternative to the tensor magnetic dif- 
fusivity employed here, one could also write down separate 
equations for the charged species (see Paper I). Assuming 
that the charged particles consist only of ions and electrons 
(denoted by the subscripts 'i' and 'c', respectively), one can 
write A = T/3i in the Hall regime, where T is the ratio of 
the Keplerian rotation time to the neutral-ion momentum 
exchange time and is the ratio of the gyrofrequency of 
charged species j to its collision frequency with the neutral 
gas (sec Appendix A). In the ambipolar diffusivity limit (not 
considered here), A = T. 

(v) e = —Vro/cs, the normalized inward radial speed at 
the mid-plane. This is a free parameter of the disc solu- 
tion, whose value is determined when it is matched to a 
self-similar global wind solution (see Section 3.5). 

(vi) ee = —VBro/cs = —cE^o/cgBz, the normalized (ver- 
tically constant) azimuthal component of the electric field, 
which measures the radial drift velocity of the poloidal mag- 
netic field lines through the disc. 



2.3 Boundary Conditions 

The complete system of equations comprises the set of ODEs 
given by equations (l)-(6), (16) and (17) as well as the 
algebraic relations (7)-(ll). This is a two-point boundary 
value problem for coupled ODEs. Eight boundary conditions 
must be formulated, either at the mid-plane or at the crit- 
ical (sonic) point of the flow, deflncd as the location where 
the vertical velocity reaches the isothermal sound speed (or 
Wz = 1). They are applied as follows. 

At the mid-plane. We begin by assigning odd symme- 
try to the radial and azimuthal components of the magnetic 
held, br and b^, which therefore vanish at 5 = 0. Consistently 

with this choice, the remaining variables (except Wz) have 
even (reflection) symmetry about the mid-plane and their 
derivatives vanish at that location. We also adopt p = 1, 
which follows directly from the normalization of p, and pre- 
scribe the radial inward velocity (measured by the parameter 
e) . The six boundary conditions applied at the mid-plane are 
thus 



bro = b^o = , 



(18) 
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po = 1 , —Wro = e ■ (20) 

Using these boundary conditions and equations (1), (2) 
and (7)-(10), we arrive at the following expressions for w^o 
and WEro'- 

2 ~2 

am e aoO"±o / \ /oi \ 

w<^o = - — ^(e-ee), (21) 

(TpO 4 ZCTpn 



the integration from the sonic point toward the mid-plane 
becomes unstable at small values of z. In carrying out this 
numerical scheme, it is essential to ensure that Wz attains a 
value of at least ~ 0.90 — 0.95 in the bisection runs before 
attempting to compute the location of the sonic point and 
the extrapolated values of the fluid variables there. If this 
is not done then the estimates of these quantities may not 
be good enough for the backward integration to the fitting 
point to be successful. 



WErO = (e-eB)H - w^o • (22) 

upo cTpo 2aQ 

At the sonic point: The sonic point is a singular point 
of equation (3), so an additional boundary condition is ob- 
tained by imposing the regularity condition 



Furthermore, the density at Zb is given by 



Ps = WzO 



(23) 



(24) 



(see the discussion just before equation 16). The density 
derivative at the sonic point can be obtained by differenti- 
ating the numerator and denominator of equation (3) (i.e. 
by applying THopitaFs rule), which yields a quadratic equa- 
tion for {dp/dz)s. Out of the two roots of this equation, we 
choose the negative one - corresponding to a positive veloc- 
ity gradient - as expected for a physical solution. 

2.4 Numerical integration of the disc equations 

The set of ODEs given by equations (l)-(6), (16) and (17) 
is integrated using the procedure first outlined in WK93. In 
brief, we start by assigning the mid-planc values of br, h^, 
p, Wr, and WEr (using equations 18 and 20-22), proceed 
to guess the value of Wzo (or, equivalently, Ps) as well as the 
position of the sonic point Zs, and then integrate from the 
mid-plane vertically upward. If the guessed Wzo is too high, 
Wz eventually diverges. If, however, the guessed value is too 
low, the vertical velocity reaches a majcimum and then be- 
gins to decrease with z. But a vertically decreasing velocity 
corresponds to a positive density gradient, which is unphys- 
ical for an isothermal disc in which the nonthermal forces 
acting in the vertical direction (namely, the vertical com- 
ponents of gravity and of the Lorentz force) tend to com- 
press the gas. This behaviour of the solutions enables us to 
bracket the correct value of Wz between two limits. We then 
improve on the guessed value using bisection until we are 
close enough to the physical solution to be able to estimate 
Za and hence (by extrapolating in z) the values of the fluid 
variables at the sonic point. The boundary conditions given 
by equations (23) and (24) can then be imposed to obtain 
the remainder of the disc variables at that location. 

To complete the derivation of the full solution, we 
renormalize z by the estimated height of the sonic point 
[z = z/zs) and integrate the equations simultaneously from 
the mid-plane (z = 0) and from the sonic point {zs = 1) to 
an intermediate fitting point (typically z ~ 0.90) while ad- 
justing the guessed values at both locations iteratively until 
the solution converges. This procedure is adopted because 



3 GLOBAL (SELF-SIMILAR) WIND MODELS 

The methodology described in Section 2 is appropriate for 

determining the initial acceleration of the wind from the 
surfaces of the disc, but it cannot be used to follow the 
evolution of the outflow on scales where the adopted thin- 
disc approximation z/r -C 1 breaks down. It is, however, 
necessary to ensure that the obtained wind solution con- 
tinues to accelerate and passes through the remaining criti- 
cal surfaces of the problem, which correspond to the Alfven 
and fast-magnetosonic critical points (e.g. BP82; Vlahakis 
et al. 2000; Ferreira & Casse 2004). Since we cannot self- 
consistently model this process in view of our radially lo- 
calized formulation of the disc structure, we approximate 
a global disc-wind solution by matching the disc solution 
to a global, radially self-similar (i.e. BP82-type) wind solu- 
tion. In our cold-wind approximation we only need to impose 
the regularity condition at the Alfven critical surface, which 
allows us to constrain one of the disc parameters (specifl- 
cally e, which quantifies the mid- plane radial velocity). We 
further simplify the treatment by assuming that ee = 0, 
which allows us to avoid relating the field-line inclination 
at the disc surface to the global magnetic flux distribution 
along the disk (see Ogilvie & Livio 2001). This is essen- 
tially the procedure employed (in the ambipolar diffusivity 
regime) by WK93. In sections 3.1-3.4 we describe, in turn, 
the governing equations, parameters, boundary conditions 
and rmmerical integration of the self-similar wind solution. 
The methodology for matching a localized disc solution to 
a global wind solution is presented in Section 3.5. 

3.1 Governing equations 

The wind is described by the steady-state, axisymmetric, 
idcal-MHD equations (e.g. Safier 1993; ?) that comprise the 
conservation of mass 



V • (pv) = 



and momentum 



pv -Vv = -VP - /aV$ -I- 



J X B 



(25) 



(26) 



for the neutral gas, the induction equation for the topology 
of the magnetic field 



Ampere's Law 



V X (u X B) = , 



J= —V X B 
47r 



(27) 



(28) 



(where, as customary, we have neglected the displacement 
current), and the solenoidal condition on the magnetic field 
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(V • B = 0). In the above expressions, P is the gas pressure 
and 3> is the gravitational potential of the central object, 

GM 



$ = 



(r2 + z2)i/2 



(29) 



where G is the gravitational constant and M is the mass of 
the protostar. In ideal-MHD flows, the magnetic field and 
velocity vectors are parallel in a frame that moves with the 
angular velocity fls of the magnetic flux surfaces, 



kB 

47rp 



+ (^B X r) . 



(30) 



where k/An is the mass load function of the wind (the ratio 
of the constant mass flux to the constant magnetic flux). 
In the poloidal (r — z) plane (subscript 'p'), this equation 

reduces to 

The variables and k satisfy (S ■ V)r2B = {B ■ V)A; = 0, 
and thus are constant along the magnetic fleld lines (or, 
equivalently, the wind flowlincs). Additional quantities that 
remain constant along the flow are the specific energy 



1 2 , ^ 



k 



(32) 



I 



- rv^ ■ 



(33) 



where h is the enthalpy per unit mass, and the speciflc an- 
gular momentum 

k ' 

which incorporates the contributions of both the matter (the 
flrst term on the right-hand side) and the magnetic fleld (the 
second term). The quantities k, e and I can be expressed in 
dimensionless form as 



and 



VGAfr^, 



(34) 



(35) 



(36) 



where the subscript 'b' denotes the location of the base of 
the wind. 

We now introduce the self-similarity scalings using the 
notation of BP82: 



z = rbX, 
Vr- = ^'(x)/(x)vKb , 

'V'l' = fl'(x)VKb , 



(37) 
(38) 
(39) 
(40) 

(41) 



In these expressions, ^' = tant^ = Br/Bz is the inclination 
of the field lines with respect to the rotation axis of the star 
and disc. At the base of the wind, we take Xb = 0, ^b = 1, 
gb = 1 and /b = 0, so the fluid velocity at the launching 
point of the outflow is exactly Keplerian. 



We now sketch the procedure followed by BP82 to ob- 
tain the set of ODEs in x that describe the self-similar wind 
solution. First, from the scalings (39)-(41), we deduce 

GM 



2 

V = 



where 



(42) 



(43) 



with the prime indicating a derivative with respect to x 
Similarly, the gravitational potential can be expressed as 

GM „ 



rb 



where the quantity S is defined by 



(44) 



(45) 



Furthermore, since we restrict our analysis to 'cold' solu- 
tions, the enthalpy term in equation (32) can be neglected 
in comparison with the other terms. Substituting equa- 
tions (32), (33), (42) and (44) into equations (35) and (36) 
yields 

2^g)-S = -l, (46) 



A: 



where we used Qb = (GM/r^y^^. The numerical value on 
the right-hand side of equation (46) is obtained by evaluat- 
ing this expression at the disk surface; it remains constant 
along the flow. 

To make further progress, one can use equation (30) to 
write 

k B4, 



9 = 



V4, 



{GM/rby/^ Anp (GM/rbY/^ 
which, together with equation (33), gives 



+ 5, 



GM ^ 



(47) 



(48) 



Then, substituting equation (48) into quation (47), one ob- 
tains 

-mA 



^(1 - m) 



(49) 



where m is the square of the Alfven Mach number (the ratio 
of the poloidal speed to the poloidal Alfven speed). Finally, 
substituting equation (49) into equation (46) and using e — 
A = —3/2, one finds 



ai - m) 



where 



T^ = C + 2S-3. 



(50) 

(51) 



Note that the gravitational plus centrifugal potential at the 

position , x} can be expressed as —{GM/rh){^'^/2 + S), 
which becomes — (3/2)GM/rb at the base of the wind. It is 
thus seen that —Tv,/2 is the gravitational plus centrifugal 
potential difference (in units of GM/vh) between the point 
, x} along a flowline and the base of the flow at = 
1 , X = 0}, so that Tw must be > for physical solutions. 

One can obtain an expression for dm/dx = x^m' 
(where x = x/xa is the vertical coordinate along a flow- 
line in units of the height of the Alfven point, subscript 
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'A') by substituting Ampere's Law (equation 28) into the 
momentum equation (26) and then combining the vertical 
component of equation (26) (with the thermal pressure term 
neglected in view of the 'cold flow' approximation) with the 
differential form of equation (50). The result, presented by 
BP82, is 



dm 
dx 

where 



= XA- 



mS'' [Bi - B2{m - 1) - Bsim - if 
^Tw(m-l)(tw-l) 



B? + BI 



Bi = 2m'x(5' - A) Jw 



B2=$(x + ^e')[(5/4)Tw + ^'-S], 

B3 = J^[xie+T^)-fix + ^0] , 



and 



= 7=-k5/ J^S — — 
-I w w 



(52) 



(53) 



(54) 
(55) 
(56) 
(57) 



(58) 



(59) 



We will also make use of 
dx 

Equation (52) has two singular points. One of them 
(m = 1) corresponds to the Alfven critical surface, which 
occurs at the location (x = Xa or x = 1) where the poloidal 
velocity component (tip) becomes equal to the poloidal com- 
ponent of the Alfven velocity (wAp). The other singular point 
{Uw = 1) occurs at the location where v-p becomes equal 
to the fast-magnetosonic speed and corresponds to the so- 
called modified fast-magnetosonic critical surface (e.g. Bo- 
govalov 1994). Although an outflow solution that is fully 
causally disconnected from the source must cross both of 
these critical surfaces (e.g. Vlahakis et al. 2000), the condi- 
tion tw is not expected to provide any additional constraint 
on our cold outflow solutions (see Ferreira & Casse 2004) 
and wo therefore ignore it in this work. 

Equations (52) and (59) can be integrated from the 
Alfven point to the base of the wind to obtain the wind 
solution. As the location of the Alfven critical point (xa) 
and the inclination of the field lines at the base of the wind 
(^b) arc not known a priori, we introduce the additional 
equations 



and 



dXA 
dx 

dx 



= 



= 0, 



(60) 



(61) 



so that these quantities can be found self-consistently when 
the equations are integrated. 

3.2 Parameters 

The global, self-similar wind solutions are specified by the 
following parameters: 



(62) 



(i) The normalized mass-to-magnetic flux ratio, 

_ 47rpUp _ ps v-Kh 

= D ~ ^ ■ 

P ^0 

(ii) The normalized total specific angular momentum, 

_ / 'rB4,\ 1 al Cs , , 

A = rvs = \- Orf,b . (63) 

(iii) The inclination of the field lines at the disc surface, 
measured by the angle that the lines make with the ver- 
tical. 



^b' = tan vPb = hrh ■ 



(64) 



The rightmost expressions in equations (62) (64) show how 
the global wind parameters can be represented in terms of 
the local disc parameters at the base of the wind (see Sec- 
tion 2.2). Note that k(A — 1) = |60b| is restricted to a limited 
range of values for physically viable solutions (see Fig. 1.2 
and WK93). For given choices of two of the above param- 
eters, the third is constrained by the requirement that the 
solution accelerates past the Alfven critical surface (see Sec- 
tion 3.5). 

3.3 Boundary Conditions 

The system of equations that describes the wind consists 
of the set of ODEs given by equations (52) and (59) (61), 
together with the algebraic equations (43), (45), (50), (51) 
and (53)-(58). This is a two-point boundary value problem; 
four boundary conditions are required, which cam be speci- 
fied either at the Alfven point or at the base of the outfiow. 
They are applied as follows. 

At the Alfven point. Two boundary conditions are im- 
posed at this location. First, by definition, 

niA = 1 • (65) 

This is a singular point of equations (50) and (52). Applying 
the regularity condition to equation (50) yields 



(66) 



Note that one can infer from equations (65) and (66) that the 

numerator of equation (52) vanishes identically at the Alfven 
point, so no additional boundary condition is obtained by 
applying the regularity condition to the latter equation. 
Applying I'Hopital's rule to equations (50) and (52), and 
using equations (65) and (66), one arrives at the following 
expressions for m'j^ and 



and 



(TwA - /1C/a)V2 
CirriA + C2m'A + C3 = , 



where 



and 



Cl = ^aT'wa(^wA — 1) , 
C2 = ^aSUxa + ^Ad)[(5/4)TwA +d- Sa] 

C3 = -4S'i^ACAXA JwA . 



(67) 

(68) 

(69) 
(70) 

(71) 
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From equations (67) and (68) one obtains {dm/dx)A = 
XAm'p^ and {d^/dx)A ~ XaCai which can be used to start 
the integration from the Alfven point toward the disc. 

At the base of the wind. Here we could, in principle, 
use the condition = 1 from the adopted self-similarity 
scaling (see equation 37). However, the integration becomes 
unstable close to the disc surface and ^' typically diverges as 
X — >^ 0. As a result, a Taylor expansion about = 1 , X = 
0} is used to obtain the fluid variables at a small distance 
{Sx w 10"'') above the disc surface. Specifically, we use 



fs. « f{,Sx + ^5x^ 



and 



where 



Jh — 



£.sx 



/b 



(3t/b-4)2 

[K2(A-l)2 + [7b]5 ' 

|(3gi+5t/b)$;/^ + I>2 
3f/b-4 



[/b 



4 



(72) 



(73) 



(74) 



(75) 



(76) 



Di = K^(A-9)(A- 1), 



(77) 

(78) 

and f/b = 1 + (.h (see BP82). The two boundary conditions 
applied at x = 5a; are then 



and 



f{6x) = fsx 



(79) 



(80) 



where the expressions on the left-hand sides of equa- 
tions (79) and (80) are obtained through the integration 

of the system of ODEs from the Alfven point toward the 
disc and those on the right-hand sides are evaluated using 
equations (72) and (73). 

3.4 Integration of the wind equations 

To integrate the wind equations, it is first necessary to 
choose the values of the free parameters k and A and sup- 
ply initial (guessed) values of xa and ^(,. In the coupled 
disc-wind solutions considered here, the adopted values of 
K and A as well as the initial estimate of ^b - are, in 
fact, calculated from the rightmost expressions in equa- 
tions (62)-(64), using the values of the fluid variables at 
the base of the wind that arc obtained at the end of the it- 
erations on the disc solution (sec Section 3.5). One can then 
use equations (67) and (68) to evaluate {d^/dx)A = xa^a 
and {dm/dx)A = XaWa and start the integration of equa- 
tions (52) and (59) from the Alfven point (a::A = 1) toward 
the base of the wind. 

Below the Alfven point (i.e. for a; < 1), the value of ^' on 
the right-hand side of equation (59) is obtained by substitut- 
ing U, 5, Tw, / and Jw (which are found from equations 43, 



45, 51, 53 and 57, respectively), as well as m and ^ (which 
are found from the integrals of equations 52 and 59) into 
equation (50). This yields the following quadratic equation 
for 

(1 - KixW)^'^ + 2xAxKi^ + (1 - Kif) 
where 



0, 



(81) 



ai-m) 



(82) 



Out of the two possible roots of this equation, we choose the 
one that satisfies the condition tan (p = ^' < tan ^ = ^/x, as 
appropriate for a solution that describes a collimating wind. 

As discussed at the end of Section 3.3, the boundary 
conditions at the disc are actually applied at a small dis- 
tance 5x above the surface, where the values derived by in- 
tegrating down from the Alfven surface are matched (via 
equations 79 and 80) to the values obtained by stepping off 
the mid-plane with a Taylor expansion. The full solution 
is then found by iterating on xa and ^f, until convergence 
is reached. Using this procedure, we derived wind solutions 
for a wide range of values of the parameters k and A (see 
Fig. 1). This solution 'library' is useful when one proceeds to 
smoothly match a radially localized disc solution to a global 
wind solution (see Section 3.5). It is available, in tabular 
form, on the VizieR data base of astronomical catalogues 
(http:/ /cdsarc. u-strasbg.fr/). 



3.5 Matching the (localized) disc and (global) 
wind solutions 

Having derived a localized disc solution as described in Sec- 
tion 2.4, one calculates the associated parameters of a self- 
similar wind solution (k, A and ^b) using the rightmost ex- 
pressions in equations (62)-(64). The parameter combina- 
tion obtained in this way does not, however, generally cor- 
respond to a wind solution that crosses the Alfven critical 
surface. In the next step, one derives a self-similar wind so- 
lution as described in Section 3.4, using the given values of 
K and A and employing the value of ^b from equation (64) 
as an initial guess in implementing the boundary conditions 
at the base of the wind. The final value of ^f, from the wind 
iteration will, in general, be different from brh, so that equa- 
tion (64) will not be satisfied. To obtain a self-consistent 
disc-wind solution, one iterates on the disc and the wind 
solutions, using e (the normalized mid-plane radial velocity) 
as an adjustable parameter, until the value of from the 
wind iteration satisfies equation (64). The CPU time to com- 
pute either a disc or a wind solution as described above is 
typically under a few seconds on a 2.4 GHz, AMD Opteron 
system. On the other hand, the CPU time associated with 
obtaining a matched disc-wind solution varies with the par- 
ticular case, but it is not unusual for the entire procedure 
to take up to 15-20 minutes. 



4 RESULTS 

In this section wo discuss the main features of radially local- 
ized solutions of well-coupled, wind-driving discs in the Hall 
diffusivity domain. In Paper I we identified four parame- 
ter sub-regimes in this domain by imposing general physical 
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Figure 1. Global, self-similar (BP82-typc) wind solutions plotted in the 5(,(= brb) ^ log-^ plane, where A is the normalized total specific 
angular momentum (including the matter and magnetic field contributions) and ^'^ measures the inclination of the field lines at the base 
of the wind. The curves are labelled by the normalized mass-to-magnetic flux ratio, k: the values of k that correspond to the darker 
curves are indicated in the figure. These solutions are available electronically in tabular form, as detailed in the text. 



constraints on viable solutions in the context of the hydro- 
static approximation, wherein the vertical velocity compo- 
nent is neglected. This approximation has made it possi- 
ble to derive algebraic relations that characterize the extent 
of each sub-regime and the distinguishing properties of the 
associated solutions. Those are summarized in Tables Al 
and A2, respectively, of Appendix A, where we also collect 
some of the algebraic expressions that were employed in the 
derivation of these results. In Section 4.1 we present repre- 
sentative solutions for these four sub-regimes (labelled by 
the Roman numerals i through iv) and compare their prop- 
erties with the predictions of Table A2. We then proceed 
(Section 4.2) to analyse the dependence of the solutions on 
the conductivity component ratio |(7h]/o-x as well as on the 
sign of the Hall conductivity cth (i.e. on the magnetic field 
polarity) and on the parameter e. We also test whether phys- 
ically viable solutions are indeed excluded from the parame- 
ter regimes that are 'forbidden' according to the hydrostatic 



analysis. Finally, in Section 4.3, wo present illustrative disc 
solutions that are matched to self-similar wind solutions and 
briefly discuss the properties of the joint disc-wind solutions 
obtained in this way. 

In all of our disc models we set the parameter eb to be 
identically zero. This parameter, which measures the radial 
drift velocity of the poloidal magnetic field lines (see Sec- 
tion 1.3.13), depends on the global distribution of along 
the disc and, therefore, cannot be obtained self-consistently 
in our localized formulation. The justification for setting 
eb = in our analysis is discussed in WK93 and in Ap- 
pendix LA. Basically, solutions characterized by the same 
value of the parameter combination (e — ee) are qualitatively 
similar, a property that results from the fact that the only 
change in the underlying system of equations brought about 
by switching to a reference frame anchored in the radially 
drifting poloidal flux surfaces (which move with the radial 
velocity vbto = — cbCs) is the substitution of Vr — vbio for 
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the radial velocity Vr in all the equations except the angular 
momentum conservation relation, which remains unchanged. 

An additional simplification that we employ is to have 
the Ohm, Pedersen and Hall conductivity terms (ao, up and 
(Jh, respectively) scale with the gas density p and magnetic 
field amplitude B as p/B^. We adopt this dependence since 
it results in the local matter-field coupling parameter A be- 
ing constant with height. This is convenient for comparison 
of our solutions with the analytic results of Paper I, which 
were effectively obtained under the same approximation, as 
well as with the solutions derived by WK93 in the ambipolar 
diffusivity regime, where the corresponding value of A was 
similarly assumed not to vary with height. Further calcula- 
tions, exploring the properties of solutions computed with 
more realistic ionization and conductivity profiles, will be 
presented elsewhere. 

4.1 Illustrative disc solutions in the Hall 
parameter sub-regimes 

In this subsection we present representative solutions for the 
four Hall parameter sub-regimes (see Table 1). We first di- 
vide them according to whether the value of the mid-plane 
Elsasser number Ao is > 1/2 (Cases i and iii; Fig. 2) or 
< 1/2 (Cases ii and iv; Fig. 3), which is one of the two clas- 
sification criteria that define these sub-regimes.^ Wc further 
divide them according to the second classification criterion, 
which is whether so = /3eoAo, the mid-plane value of the 
ion slip factor, is > 1 (Cases i and ii, shown in the left pan- 
els of Figs. 2 and 3, respectively) or < 1 (Cases iii and iv, 
depicted in the corresponding right panels of these figures). 
The variables /3io and /Seo are, respectively, the mid-plane 
values of the ion and electron Hall parameters; see equa- 
tion Al. Now, the hydrostatic analysis implies that viable 
solutions must satisfy so > 1 in the ambipolar diffusivity 
regime (WK93) and Sq < 1 in the Ohm diffusivity regime 
(Paper I). One can therefore classify Cases (i) and (ii) as 
being in the 'ambipolar diffusion-modified' Hall regime and 
Cases (iii) and (iv) as being in the 'Ohm diffusion-modified' 
Hall regime. This description is supported by the fact that 
the parameter constraints and solution properties derived 
for the Hall Case (i) are identical to those in the ambipolar 
diffusivity limit, with a similar correspondence holding be- 
tween the Hall sub-regime (iii) and the Ohm sub-regime (i) 
(see Paper I). 

We use the solution for Case (i) (left panel of Fig. 2) 
as a representative model that illustrates the overall prop- 
erties of well-coupled wind-driving discs. The top panel of 
this figure shows the gas density as well as the radial and 
azimuthal components of the magnetic field as functions of 
height between the midplane and the sonic point. The bot- 
tom panel depicts the velocity components. This layout is 
used also in all the other figures in this paper that depict 
the vertical structure of the disc solutions. 

The representative disc solution exhibits the following 
overall properties (WK93). In the quasi-hydrostatic layer 
straddling the mid-plane (0 < \z\ < 0.4 in this example), 

^ Computing solutions for Ao < 1/2 has proven to be numerically 
challenging. For this reason, the solutions presented in Fig. 3 cor- 
respond to values of Aq that are only slightly smaller than 1/2. 



the magnetic field lines are radially bent and azimuthally 
sheared (see equations 5 and 6), and matter loses angular 
momentum to the magnetic field. This enables the disc mate- 
rial to flow toward the protostar (i.e. e s —Wro > 0) and re- 
sults in a radial drag on the ionized disc component (assum- 
ing that cb remains < e; see WK93 and Appendix LA). This 
drag must be balanced by magnetic tension (see equations 5, 
8 and 10), which is consistent with the field lines bending 
radially outward. Furthermore, as the magnetic tension also 
partially supports the fluid against the gravitational pull of 
the protostar (see equations 1 and 5), the g£is motion in 
this region is sub-Keplerian. Further up but still below the 
base of the wind (defined as the height Zb whore — vk 
and located in this solution at z ~ 3.0) lies the transition 
zone (0.4 < |z| < 3.0). In this layer the magnetic pres- 
sure comes to dominante the thermal pressure because of 
the strong drop in the gas density away from the midplane, 
and the magnetic field lines are locally straight. Note that 
everywhere within the disc the radial velocity is negative, 
the azimuthal velocity is sub-Keplerian, and the vertical ve- 
locity is subsonic (highly so in the quasi- hydrostatic layer). 
Finally, as the angular velocity of the magnetic field lines 
is nearly constant along the poloidal flux surfaces (it is ex- 
actly constant only when ee is identically zero; see Konigl 
1989), their azimuthal velocity increases with radius until, 
eventually, they overtake the fluid (whose azimuthal veloc- 
ity decreases with radius according to the Kcplcrian rotation 
law). In the hydrostatic limit (vz — > 0) and with eb = this 
occurs at the height Zb, and under the same assumptions 
this is also the location where Vr changes sign.^ The solu- 
tion region above Zb (between Zb w 3.0 and the sonic surface 
at Zs ~ 4.7) represents the base of the wind. In this region Vr 
is > 0, u^f, is supcr-Kcplcrian and the magnetic field transfers 
angular momentum back to the matter, causing the upper- 
most layers of the disc to be driven out centrifugally (with 
the mass flux in the wind effectively fixed by the regularity 
condition at the sonic critical point of the disc solution). 

The distinguishing properties of the four Hall sub- 
regimes, as inferred in the context of the hydrostatic ap- 
proximation, are listed in Table A2. To see how the expres- 
sions derived in this approximation for \dbr/db^\a, h = h/hT 
(where h is the magnetically reduced density scaleheight, de- 
fined as the location where the density drops to p^j \fe) and 
Zbjh compare with the properties of the full solutions listed 
in Table 1, we need to relate the parameters /3io, /3eo and To 
used in the two-component plasma formulation of Paper I to 
the tensor conductivity components employed in this work. 
The parameter To is the mid-plane value of the neutral- 
ion coupling function given by equation A2 and is equal to 
the Elsasser number in the ambipolar diffusivity limit. The 
relationships between these variables and the conductivity 
ratios can be obtained from equations 1. 95-1. 97 and reduce, 
in the limit q = Pi/ <C 1, to 

To«Aofl^^-?iV' , (83) 



In an e-xajct solution, where ii^ is > near the top of the disc, 
the heights where = v^i, = rCl^ and u,. = do not exactly 
coincide, although typically they remain close to each other. 
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Figure 2. Structure of illustrative disc solutions as a function of height above the midplane (measured in units of the tidal scaleheight 
Ht) for the Hall conductivity sub-regimes labelled in Table Al as Case (i) (left panels) and Case (iii) (right panels) and the model 
parameters shown in the figure (where all the listed conductivity values pertain to the mid-plane). For these two cases, the matter-field 
coupling parameter Aq > 1/2. In both solutions, the top panel shows the gas density (normalized by the midplane value) as well as 

the radial and azimuthal components of the magnetic field (normalized by the vertical component, which is constant with height). The 
bottom panel depicts the velocity components in a frame rotating at the Keplerian speed, normalized by the isothermal sound speed Cs. 
The curves terminate at the sonic point. The key properties of both solutions are listed in Table 1. 
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Figure 3. Same as Fig. 2 except that the depicted solutions correspond to Cases (ii) (left panels) and (iv) (right panels) of Table Al, 
for which Aq < 1/2. Key properties of these solutions are summarized in Table 1. 
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Table 1. Key properties of the solutions shown in Figs. 2 and 3. We list the lst-3rd inequalities of Table Al (representing, respectively, 
the requirement for a sub-Keplerian flow within the disc, the wind-launching condition and the lower bound on the height of the base 
of the wind) for each case as well as the mid-planc values of the ion slip factor s = /3e/3i and of the matter— field coupling parameter 
A = CTx, which is equal to T|/3i| in the limit (|(ThI/5'± — 1) "C 1 (see equation 1.95). These parameters form the basis of the classification 

scheme of viable solutions in the Hall diffusivity regime (sec Section 1.5). The base of the wind and the sonic-surface heights (z\y and 
2a, respectively) arc in units of both the tidal scalcheight h'Y and the actual (magnetically reduced) pressure scaleheight h, where the 
latter is determined in each case as the value of z where the density drops to po / y/e. The mass accretion rate per disc circumference is 
evaluated from Mi^/2iTro = ~2 J^^^ pVrdz, and the listed value (A4in) is this quantity normalized by poCg/ix. 

Illustrative Case (i) Case (iii) Case (ii) Case (iv) 

disc solution Fig. 2 (left) Fig. 2 (right) Fig. 3 (left) Fig. 3 (right) 
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V2To ~ " ~ ~ 
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0.2 < 1 < 2 < 9 


0.1 < 1 < 2 < 10 


0.4 < 0.4 < 1 < 3 


3 < 4 < 1 < 3 
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and 
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5-± 



(84) 



(85) 



4.2 Parameter dependence of the solutions 

4-2. 1 Dependence on the ratio \aYi\/d± 

The conductivity ratio crH/5-± serves to distinguish the am- 
bipolar diffusivity domain (ao 3> 5"p 3> j(5"H|) from the Hall 
diffusivity regime (iTp <^ \dYi\ <^ 5"o): the ambipolar diffu- 
sivity limit corresponds to this ratio tending to 0, whereas 
the Hall limit corresponds to this ratio tending to 1. To il- 
lustrate the difference between these two regimes, it may be 
instructive to compare a solution in the ambipolar diffusivity 
regime, which satisfies Aq > 1 and sq > 1 (see Paper I), with 
a solution in the Hall sub-regime (i), which, as we already 
noted in Section 4.1, satisfies the same inequalities and, in 
fact, identical parameter constraints (first line in Table Al). 
We show such a comparison in Figs. 4 and 5. Fig. 4 depicts 
the vertical run of the density as well as of the transverse 
magnetic field and the velocity components for two represen- 



tative solutions in the above-mentioned diffusivity regimes 

for the same values of the parameters Ao, ao and e. Fig. 5, in 
turn, shows the dependence of the density decrement ps/po 
at the sonic surface and of the characteristic heights Zh = h, 
5b and Zs on the parameter e for the same two choices of the 
conductivity components (as well as for a solution with an 
intermediate value of crn/ax = 5/-\/2) and the same values 
of Ao and ao as in Fig. 4. 

The most noticeable differences between the two solu- 
tions shown in Fig. 4 are the faster change of the flow vari- 
ables (corresponding to a smaller value of h) and the larger 
vertical extent of the displayed solution (corresponding to 
a higher value of Zs) in the Hall case. The higher value of 
Zs implies a lower density at the sonic surface and there- 
fore a lower normalized upward mass flux (= ps) for the 
Hall solution (see Fig. 5). The decrease (increase) of Zh (zs) 
with increasing |o"h|/(3"± when all the other parameter val- 
ues are held constant is also evident from the curves plot- 
ted in Fig. 5. The latter figure further demonstrates that 
the behaviour of Zs is similar to that of 2b, which suggests 
that Zb could be used as a proxy for Zs in analysing this 
trend. The parameter dependence of h and Zb can be es- 
timated using the hydrostatic approximation, as discussed 
in Paper I, and is given by equations (A5) and (A6), re- 
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Figure 4. Structure of a wind-driving disc solution as a func- 
tion of 2/ /ix for the disc parameters shown in the figure and two 
different configurations of the conductivity tensor (whose listed 
components all pertain to the mid-plane). The thick lines corre- 
spond to the ambipolar diffusivity limit (for which cth = and 
ctq S> CTp) and the thin lines to a case where cth > o^p. The 
normalized upward mass flux is pwz = Ps = 1-1 x 10~^ for the 
solution in the ambipolar diffusivity limit and 1.4 x 10~^ in the 
Hall case. 



spectively. Given that Ao To|Ao| as IcthI/^x — >■ 1 and 
assuming that Ao, ao and e are held constant, the above 
equations imply that, to leading order in |/3io|, h oc |/?io| and 
Zh oc in the ambipolar regime and in the Hall sub- 

regimes (i) and (ii), for which so = /SeoAo ^ 1- On the other 
hand, equation (84) impUes, given that a± = Ao /oq (equa- 
tion 1.87) and ao are fixed, that |/3io| decreases as |i5-h|/5"± 
goes up (and crp/ax = (o"5^ — Sr'^Y^^ /a±_ goes down). Taken 
together, these results explain the manifested dependence of 
the solutions in Figs. 4 and 5 on IcthI/ctx. 

As a further check on the apphcabihty of the analyt- 
ical framework developed in Paper I, we note that equa- 
tions (A5) and (A6) imply a similar dependence of 2h and 
Zh on |aH|/(5"x in the limit so ^ 1 (i.e. in the Ohm diffusion- 
modified Hall regime; Cases iii and iv). In fact, fixing the 
values of the same parameters as before, these two equations 
imply that, to leading order in small ratios, h oc l/|/3o()l and 
Zb oc |/9eo| in this case, while equation (85) shows that |/3eo| 
increases as |5'h|/5'x goes up (assuming again that d±_ and 
do remain unchanged). We have verified that solutions in 
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Figure 5. Top panel: Density at the sonic point {ps = Ps/po) as 
a function of the normalized inward velocity at the midplane (e) 
for three configurations of the conductivity tensor characterized 
by the same values of ctq and a± but different values of the ratio 
(7li/a±. (All the listed conductivities pertain to the mid-plane.) 
The other model parameters have the same values (indicated in 
the figure) in each case. Bottom panel: As above, but now plotting 
the vertical location of the sonic point (is), the base of the wind 
(ib) and the magnetically reduced scaleheight (zh = h). In both 
panels, the solid lines show the ambipolar diffusivity limit (cth = 
0), the dashed lines represent the case where cth = crp and the 
dot-dashed lines depict a case where cth > crp. 



these regimes indeed exhibit the expected dependence on 

\d-u\/d-±. 

Yet another test of the predictions of the hydro- 
static analysis regarding the dependence of the solutions on 
|5"h|/5"x can be constructed using equation (83) for the pa- 
rameter To. As discussed in Section 1.6 and illustrated in 
Fig. 1.2, this analysis indicates that the requirement To ^ 1 
(i.e. that the mid-plane neutral-ion momentum exchange 
time be shorter than the local orbital time) is a funda- 
mental constraint on viable wind-driving disc models of the 
type that we consider, and applies to solutions in all the 
diffusivity regimes (ambipolar. Hall and Ohm). This con- 
dition follows directly from the inequalities To|/3io| > 1/2 
and |/3io| < 1 that characterize the Hall Cases (i) and (iii), 
but, as can be seen by combining the first two parameter 
constraints reproduced in the first row of Table 1, it can be 
formally inferred (within the framework of the hydrostatic 
analysis) to apply also in the other Hall sub-regimes. The 
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Figure 6. Dependence of the ratio cth/ctx on the parameter To 
for solutions corresponding to the indicated parameter values. 
The solid and short-dashed curves show solutions for the Hall sub- 
regimes (ii) (eomputcd with ctq = 41 at the mid-plane) and (iv) 
(computed witii oq = 20 at the mid-plane), respectively. In both 
cases, (Txo = 3.1. All viable solutions lie to the right of the vertical 
long-dashed line, where To > 1, as predicted by the hydrostatic 
analysis. 

illustrative solutions for the Cases (ii) and (iv) in Fig. 3 
satisfy this inequality, as verified explicitly in Table 1. One 
can further test the foregoing result in these sub-regimes by 
holding Ac, ax and bo fixed (as was done above) and de- 
creasing IcthI/ctx (or, equivalently, increasing ap/d_L) until 
(if Aq is sufficiently small) To declines to a value < 1 where, 
according to the hydrostatic analysis, viable solutions cease 
to exist. Such calculations are presented in Fig. 6, where we 
plot lanl/o-x as a function of To for representative solutions 
in these two sub-regimes. It is seen that all viable solutions 
are indeed restricted to the region Tq > 1, validating the 
above prediction. 

4.-2.2 Dependence on the field polarity 

As was noted in Section 1.2.1, the Ohm and Pedersen con- 
ductivities arc always positive, even under a global reversal 
of the field polarity, since gq is not a function of the mag- 
netic field strength and ap only contains magnetic terms 
that scale as . The Hall conductivity, however, has an 
overall linear dependence on B = |B| sgn{Bz} and can thus 
assume both positive and negative values depending on the 
direction of the vertical field component. The dependence 
of (Th on the magnetic field polarity was shown in Paper I 
to affect both the extent of the parameter ranges where vi- 
able wind-driving disc solutions can exist in the Hall domain 
and the properties of these solutions. We now briefly sum- 
marize these results, which were obtained in the hydrostatic 
approximation assuming an ion-electron plasma. 

The range of values of the parameter f5 = l//3io for 
viable solutions was found to be restricted by the following 
two conditions (see Section 1.6): 

(i) Sub-Keplerian flow below the base of the wind. This 
implies that the normalized azimuthal velocity Wtf, is < 
within the disc, which translates into d6r-/d6,^ < and /? > 
-2To (equations I.106-I.108). 




Figure 7. Normalized sonic-point density ps as a function of the 
parameter Tq for solutions in the Hall sub-regime (i), illustrating 
the dependence of viable disc models on the magnetic field polar- 
ity (i.e. the sign of cth)- Curves for cth > are shown by light lines 
and for <th < by heavy lines, and are labelled by the normalized 
mid-plane inflow velocity e. (All the listed conductivities pertain 
to the mid-plane.) The value of |/3| s l/|/3io| is the same (1.45) 
for all the solutions. The results are broadly consistent with the 
prediction of the hydrostatic analysis (equation 87) that there 
should be no viable solutions for an < when To drops below 
2|/3| (i.e. to the left of the vertical dashed line). 

(ii) Sub-Keplerian azimuthal velocity of the magnetic flux 
surfaces, or lOBro = —r{^Bo — Q,iC)/ca > 0, where fieo = 
—cEro/rBo is the angular velocity of the flux surfaces at 
the disc mid-plane. This implies (see equation 1.113) that 
either P < — 2To or /3 > —To/2 in the Hall domain. 

The above two constraints together imply 

/? > — To/2 in all the Hall sub-regimes. (86) 

By combining this inequality with the classification criterion 

T()/|,5| > 1/2 that distinguishes the Hall sub-regimes (i) 
and (iii) (see Section 4.1), one obtains 

-To/2 < /? < 2To Cases (i) and (iii) . (87) 

Equation (87) predicts that, although there could be 
both positive- and negative-polarity solutions in the Hall 
sub-regimes (i) and (iii), no viable solutions should exist 
in these cases when /3 decreases below — To/2. This pre- 
diction of the hydrostatic analysis is examined in Fig. 7, 
which plots Ps, the normalized density at the sonic point, as 
a function of the coupling parameter To for Case-(i) solu- 
tions derived using the indicated model parameters and cor- 
responding to both positive and negative values of /3 oc cth 
(sec equation 84). Since the magnitudes of all the conduc- 
tivity tensor components are held constant, all the solutions 
are characterized by the same absolute value of the ion Hall 
parameter (]/3io| ~ 1/1.45). The figure verifies that viable 
solutions cease to exist in the region to the left of the ver- 
tical dashed line, which corresponds to |/3| > To/2, when 
P < 0.^ 

^ Note that each of the /9 < solution curves in Fig. 7 termi- 
nates at a finite value of Tq that is lower the higher the value 
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Figure 8. Comparison of two disc solutions that differ only in the 
sign of the Hall conductivity; (th > for the solution shown by 
thin lines, and o-jj < for the solution depicted by thick lines. The 
normalized upward mass flux pWz is 9 X 10~^ in the (Th > case 
and 5 x 10~^ in the cth < solution. All the listed conductivities 
pertain to the mid-plane. 



By similarly combining the inequality (86) with the 
classification criterion To/|/3| < 1/2 that distinguishes the 
Hall sub-regimes (ii) and (iv), one obtains 



/? > 2To Cases (ii) and (iv) , 



(88) 



which implies that, for these sub-regimes, self-consistent so- 
lutions exist only if the field has a positive polarity. We 
verified this prediction by adopting the values of ao , e and 
a± specified for Cases (ii) and (iv) of Tabic 1 and chang- 
ing the value of o-h/ctx: we found that no negative-polarity 
solutions could be obtained in these cases. 



of e. To understand this behaviour, note that, given the param- 
eters that are held fixed in this figure, To scales as Oq. For a 
fixed value of ao, ps increases as e goes down (as seen also in the 
top panel of Fig. 5). There is a maximum sonic-point density ps 
that can be attained for the chosen value of ay, corresponding to 
the maximum outflow rate for a consistent solution, and this, in 
turn, determines the magnitude of e for the solution curve that 
terminates at the given value of ag. A lower value of ao corre- 
sponds to a higher lower bound on e because it implies a higher 
value of |6,^bl (see equations A3, A4 and 83) and hence a stronger 
torque on the disc and a correspondingly higher mass accretion 
rate (reflected in the value of e). 
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Figure 9. Variation of the normalized sonic-point density (upper 
panel) and of the chara<;teristic heights Zb , % z-^ = h (lower 
panel; cf. Fig. 5) with the model parameter e for the same con- 
ductivity tensor components employed in Fig. 8. The solid lines 
correspond to (jh > and the dashed lines to (jh < 0. 



Even when both positive and negative values of (3 are 
allowed, as is the case in sub-regimes (i) and (iii), changing 
the field polarity modifies the properties of the solutions. 
This point, which was already made in WK93, is illustrated 
in Figs. 8 and 9, which compare disc solutions that differ 
only in the sign of an.'' The main difference between the 
two solutions in Fig. 8 is evidently the value of the sonic- 
point height Za (corresponding to the termination point of 
the depicted solution), which is larger (~ 4.3) in the positive- 
polarity case than in the negative-polarity solution (~ 2.4). 
Fig. 9 confirms this trend: The wind-launching surface (ib) 
and the sonic surface (zb) are located higher above the mid- 
plane, and the mass outflow rate (measured by ps) is cor- 
respondingly lower, in the positive-polarity solutions. The 
value of £h = h, which represents the scale on which the 
density and transverse magnetic field components undergo 
their strongest variation, is also different between the two 
sets of solutions (it is larger for the negative-polarity model), 



In discussing the differences between Hall-domain solutions 
with opposite field polarities in Section 1.6 we explicitly referred 
to sub-regime (iii), whereas the solutions shown in Figs. 8 and 9 
correspond to sub-regime (i). It should, however, be clear from 
the present discussion that the behaviour of solutions in these 
two sub-regimes is similar in this regard. 
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although this difTerence is modest in comparison with the 
change in Zh and Zs, and it goes in the opposite direction. 

The above behaviour can be understood using the ana- 
lytic expressions derived in paper I, and reproduced in Ap- 
pendix A. Employing the same reasoning as in our analysis 
of Figs. 4 and 5 in Section 4.2.1, we employ the height Zh 
of the base of the wind as a proxy for Zg in this analysis. 
Using equations (A5) and (A6), we estimate that, to lead- 
ing order, h+/h- ^ (2To - l/?!)/(2T„ + and 2b+/2b- « 
{Tl + 5T„|/?|/2 + l3'^)/(rl- 5To|/?|/2 + 13''), whore the sub- 
scripts and '-' refer to the positive- and negative-polarity 
cases, respectively. Using also equations (83) and (84), we 
infer that, for the parameters adopted in Figs. 8 and 9, 
h+/h- ~ 0.85 and ib+Z^b- ~ 2.28. These estimates are 
entirely consistent with the numerical results (for a fixed 
value of e) shown in the two figures. Furthermore, from equa- 
tion (A3) we confirm that |6rb/6<Ab| is 1 and hence (us- 
ing equation A4) that b^h ~ V^/cio, independently of the 
sign of /3, reproducing the actual behaviour of the solutions 
in Fig. 8. Setting &rb+ ~ Kh-, we then infer from equa- 
tions (A3) and (A5) that bri,h+/b^h- ~ h+/h- (« 0.85 for 
the adopted parameters). This conclusion, too, is consistent 
with the result exhibited in Fig. 8. 

Although the specific angular momentum at the base 
of the fiow is mostly magnetic (corresponding to the wind 
model parameter A being ^ 1), the fact that 6</,b+/&<^b- < 1 
does not imply that the value of A is larger in the negative- 
polarity case. In fact, the converse is true, as can be seen 
using equation (63), from which it follows that A+/A- w 
(6,^b+/b0b-)(Ps-/Ps+)- Given that ps-/Ps+ is typically » 
b^b- /b^b+ (as the solutions in Fig. 8 demonstrate), wo find 
that A+/A_ 3> 1. Physically, the magnetic torque acting on 
the disc (oc B-^B^) and correspondingly the mass accretion 
rate and the rate of inward angular momentum advection 
are slightly larger in the negative-polarity case. However, 
the mass outflow rate is significantly larger in this case and 
therefore the specific angular momentum (the angular mo- 
mentum per unit mass) has to be rrmch smaller in order for 
the rate of inward (radial) and outward (vertical) angular 
momentum transport to balance each other. An alternative 
way of arriving at this conclusion is to consider the wind 
solutions that are self-consistently matched to these disc so- 
lutions. The matched wind solutions lie on a ^b' = const 
curve in the « — A wind parameter space (see equation 64) , 
and along such a curve higher values of /t oc Ps (see equa- 
tion 62) correspond to lower values of A (see Fig. 2 in BP82). 

4-2.3 Dependence on the radial velocity parameter e 

To illustrate the dependence on the parameter e (the normal- 
ized mid-plane inflow velocity), we plot in Fig. 10 two solu- 
tions in the Hall sub-regime (i) that differ only by the value 
of this parameter. It is seen that, as e decreases from 1.0 
(thick linos) to 0.5 (thin lines), the scaloheight h (the scale 
on which the density and transverse magnetic field compo- 
nents vary most strongly) increases and Zs (the vertical ex- 
tent of the displayed solution, corresponding to the location 
of the sonic surface) decreases. These trends are also evident 
in the solutions depicted in Figs. 5 and 9. As we did in Sec- 
tions 4.2.1 and 4.2.2, we use Zb as a proxy for Zs and employ 
the hydrostatic-approximation equations (A5) and (A6) to 
analyze this behaviour. These two expressions show that. 
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Figure 10. Vertical structure of a wind-driving disc solution for 

two different values of the radial velocity parameter: e = 1 (thick 
lines) and e = 0.5 (thin lines). The other model parameters are 
listed in the figure (with all the conductivity values pertaining to 
the mid-plane). 



with all the other parameters remaining unchanged, h oc 1/e 
and Zb oc e, which can directly account for the exhibited 
trends. In both cases, the dependence on e can bo traced 
to the scaling oc Wr implied by the angular momentum 
conservation relation (equation 2) in the hydrostatic limit 
(see Section 1.4.3). 

The hydrostatic analysis presented in Paper I also leads 
to a lower limit on the value of e, derived from the require- 
ment that the base of the wind be located above a density 
scaleheight, i.e. Zb > h. This implies eTo > (18)'^'''* « 2 
(see Table A2) and represents the third constraint listed in 
Table Al. To verify the applicability of this condition, we 
constructed a solution using the same parameters as in the 
illustrative Case-(i) solution depicted in Fig. 2 except that 
we chose a smaller value of e so that the above constraint 
is no longer satisfied. The result, shown in Fig. 11, confirms 
that, when the above inequality is violated, the solution is 
no longer physically viable. This is evidenced by the down- 
ward turn of the |6^| and br curves above a certain height. 
A similar solution was presented in Fig. 6 of WK93, who 
explained the origin of this behaviour by noting that, when 
e decreases to a sufficiently low value, Zs becomes so small 
and (correspondingly) ps so large that the upward mass flux 
carries more angular momentum than is brought in by the 
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Figure 11. Same as the left-hand side of Fig. 2, but for e = 0.2. 

In this case To = 9.1 (sec Table 1), so the inequality eTg > 2 (the 
third constraint on viable solutions in the Hall sub-regime (i); see 
the first row in Table 1) is not satisfied. As discussed in the text, 
the resulting disc solution is unphysical. 



accretion flow. Consequently, the gradient of changes sign 
(with that of br following suit) as the magnetic field starts 
depositing angular momentum back into the flow even before 
the nominal top of the disc is reached. Such a configuration 
is likely unstable (e.g. Cao & Spruit 1994). 



4.3 Matched disc— wind solutions 

The solutions presented in Sections 4.1 and 4.2 have involved 
the disc model alone. This was done in order to simplify the 

derivations and was an adequate approach given that we 
were primarily interested in studying the constraints on the 
disc model parameters. However, to justify this treatment, 
it is necessary to demonstrate that matched disc-wind so- 
lutions can be obtained for similar disc parameter combi- 
nations. This is done in this subsection, where we present 
two examples of self-consistently matched (local) disc and 
(global) wind solutions. These solutions correspond to the 
same values of the parameters Ao and oo and of the abso- 
lute values of the conductivity tensor components but have 
opposite signs of an- they were obtained using the proce- 
dure described in Section 3.5 and arc shown in Figs. 12 
(the positive-polarity case) and 13 (the negative-polarity 
case). As is evident from equations (62) and (63), combin- 
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Figure 12. Illustrative disc-wind solution obtained by matching 
a radially localized disc solution to a global (self-similar) wind so- 
lution using the procedure outlined in Section 3.5. The disc model 
parameters are shown in the figure (with all the conductivity val- 
ues pertaining to the mid-plane) whereas those of the matched 
wind solution are /t = 7.9 x 10"*, A = 174 and = 1.42. This 
solution satisfies the constraints specified in Table Al for the Hall 
conductivity sub-regime (iii): 0.1 < 0.97 < 2 < 6.4 < 10. 



ing the two solutions requires a specification of the param- 
eter Cs/wk ~ h-rlr: the above examples were computed for 
c/wK = 0.05. 

Our matched solutions correspond to the Hall sub- 
regime (iii) and are seen to be very similar to the analogous 
solutions we obtained by considering the disc alone (right 
panels of Fig. 2). The qualitative differences between the 
cth > and cth < solutions (essentially a higher value of 
2s and lower values of h and of in the positive-polarity 
case) are also the same as those found in the 'windless' solu- 
tion (Figs. 8 and 9). This is as expected, since the matching 
to the wind solution does not affect the physical require- 
ments on a viable disc solution; rather, what it does is fix - 
through the imposition of the Alfven regularity condition on 
the wind solution - the value of one of the disc model param- 
eters that is chosen arbitrarily when the wind solution is not 
taken into account. In our treatment (see Section 3.5), the 
disc parameter that is determined in this way is e: it is equal 
to 0.68 in the cth > solution shown in Fig. 12 and to 0.98 
in the ctr < solution shown in Fig. 13. In an even more 
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Figure 13. Same as Fig. 12, but with cth having the opposite 
sign. All the other disc model parameters except for e, which is 

determined by imposing the Alfvcn regularity condition on the 
wind solution, are the same as in Fig. 12. In this case the param- 
eters of the matched wind solution are k = 1.2 X 10^'^, A = 144 
and = 1.43. This solution, too, satisfies the hydrostatic- 
analysis constraints for the Hall conductivity sub-regime (iii): 
0.1 < 0.97 < 2 < 9.2 < 10. 



self-consistent solution, the value of the parameter ee would 
also be determined by taking into account the conditions 
outside the disc - in this case the magnetic flux distribution 
along the disc surface (see Section 1.3.1.3) - rather than ar- 
bitrarily setting it equal to zero as has been done in this 
paper and in the analytic derivations of Paper I." 



5 CONCLUSION 

In this paper we continue our study of the viability and 
properties of weakly ionized protostellar accretion discs that 
transport their excess angular momentum vertically through 
the surfaces by means of centrifugally driven winds. In view 
of the suggestive evidence that this situation is realized in 



^ A self-consistent treatment along these lines has, in fact, al- 
ready been implemented in the fully global self-similar disc/wind 



model constructed by Teitler (2010). 



at least some protostellar systems (e.g. Ray et al. 2007), 
and yet mindful of the fact that the total radial extent of 
such wind-driving disc regions is still unknown, we consider 
radially localized disc configurations (which, however, are 
joined to a global wind model). We employ the formulation 
devised by WK93 for modelling steady, geometrically thin, 
vertically isothermal and nearly Keplerian discs in which 
magnetic diflfusivity counters the shearing and advection of 
the magnetic field. WK93 assumed that the charge carri- 
ers were singly-charged ions and electrons and focussed on 
the ambipolar diffusivity regime. In Paper I we extended 
this model to the Hall and Ohm diffusivity regimes using 
the conductivity-tensor formalism, although we reverted to 
using the multifluid approach employed by WK93 for deriv- 
ing parameter constraints on physically viable solutions in 
these regimes. This derivation generalized the correspond- 
ing results of WK93 for the ambipolar diffusivity regime 
and was similarly carried out in the context of the hydro- 
static approximation, in which the vertical velocity compo- 
nent is neglected inside the disc (which results in several of 
the differential equations for the disc structure simplifying 
to algebraic relations). 

The hydrostatic analysis of paper I indicated that vi- 
able wind-driving disc solutions correspond to four param- 
eter sub-regimes in the Hall diffusivity domain and three 
in the Ohm domain. It also led to analytic estimates of 
the magnetically reduced (due to magnetic pressure-gradient 
'squeezing') density scale height and of the location of the 
disc's surface (where the inflow turns into an outflow), £is 
well as of other pertinent quantities. These results are sum- 
marized in Appendix A. In this paper we test these predic- 
tions by constructing exact solutions of the disc equations. 
We concentrate on the Hall regime in view of its expected 
importance in the inner regions of real systems; we do not 
consider the low-ionization Ohm regime in this paper given 
that it may have limited relevance to wind-driving proto- 
stellar discs. We characterize the solutions in terms of the 
conductivity-tensor components (i.e. the Pedersen, Hall and 
Ohm conductivities). However, to facilitate the comparison 
with the analytic results of Paper I, which were derived in 
the framework of the multifluid formulation, we assume that 
the ratios of these terms arc constant with height in the disc 
and, more speciflcally, that they scale with the density and 
field amplitude as p/B^, which implies that the matter-field 
coupling parameter (the Elsasser number A) is also constant 
with height. Wo require the derived solutions to cross the 
sonic critical surface but we do not contirmo the integration 
past that surface; this is sufficient for the comparison with 
the analytic results and greatly simplifies the calculations. 
However, as recapitulated below, we also demonstrate that 
these solutions can be matched to wind solutions that ex- 
tend to large distances (and, in particular, cross the Alfven 
critical surface). Our findings can be summarized as follows. 

(i) Our numerical solutions are in broad agreement with 
the parameter constraints obtained under the hydrostatic 
approximation for the four Hall parameter sub-regimes. In 
the regions of parameter space that are excluded by the 
above constraints, wind-driving disc solutions cannot be ob- 
tained or are unphysical. 

(ii) All viable solutions satisfy the constraint To > 1 (see 
equation A2). Physically, this condition expresses the re- 
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quiremcnt that the mid-plane neutral-ion momentum ex- 
change time be shorter than the disc orbital time. As dis- 
cussed in Paper I, this requirement is predicted to apply 
in each of the diffusivity regimes and is evidently a funda- 
mental constraint on disc solutions of the type considered 
here. 

(iii) For the same values of ao, Aq and e (the normalized 
mid-plane magnetic field amplitude, matter-field coupling 
strength and radial speed, respectively), increasing the rela- 
tive contribution of the Hall conductivity (the ratio |(Th1/5"±) 
results in a smaller (magnetically reduced) density scale- 
height h) and in larger heights of the disc and sonic surfaces 
(2b and Zs, respectively). A higher value of Js in turn implies 
a lower density at the sonic surface (ps) and hence a lower 
mass outflow rate. 

(iv) The magnetic field polarity affects both the proper- 
ties of the solutions and the extent of the parameter regimes 
where viable solutions exist when the Hall current is dynam- 
ically important, reflecting the dependence of the Hall con- 
ductivity on the sign of Bz- Specifically, we confirmed the 
following dependence of the solutions on the field polarity 
(with a positive polarity corresponding to being parallel 
to the disc rotation vector): 

(a) Hall sub-regimes (i) and (ni). The parameter /3, 
which is equal to the inverse of the mid-plane ion Hall 
parameter (equation Al) and thus scales inversely with 
the signed magnetic field amplitude B, is predicted to lie 
in the range — To/2 < l3 < 2To. We verified that no solu- 
tions exist for /? < — To/2 and that positive- and negative- 
polarity solutions have distinct properties. In particular, 
when the sign of an is changed from < to > and all 
the other parameter values remain the same, 2b and 2s 
are increased and the wind outflow rate correspondingly 
goes down. 

(b) Hall sub-regimes (ii) and (iv). These parameter 
regimes are characterized by ,8 > 2To. We verified that 
only positive-polarity solutions can be obtained in this 

case. 

(v) Decreasing the inward radial speed (e) increases the 
scale over which the fluid variables vary most strongly (the 
magnetically reduced density scaleheight h) but decreases 
both 2b and 2s. Furthermore, when the lower limit on the 
value of e, obtained from the requirement that 2b > ^ (the 
third inequality in Table Al) is violated, the computed so- 
lutions are found to be unphysical (in that hr and \b^\ start 
to decrease with 2 above a certain height). 

We also detail the procedure for obtaining global (radi- 
ally self-similar) 'cold' wind solutions following the method- 
ology introduced by BP82. Wc compute solutions of this 
type for a large range of values of the wind model parame- 
ters K, A and ^f, (the normalized mass-to-flux ratio, specific 
angular momentum and fleld-line inclination at the disc's 
surface, respectively). Tables of these solutions are avail- 
able on the VizieR data base of astronomical catalogues 
(http://cdsarc.u-strasbg.fr/). As our radially localized and 
geometrically thin model cannot be used to follow the prop- 
agation of the outflow far from the disc, wc match our disc 
solution to a BP82-type wind solution by adjusting one of 
the disc model parameters (e) and iterating on the disc 
and wind calculations until the full solution converges. We 



present illustrative solutions of this type that demonstrate 
that matched disk/wind configurations can be obtained for 
parameter values that are very similar to those of the merely 
transonic solutions employed in our parameter-space analy- 
sis. 

The accretion process in protostellar discs may in- 
volve a variety of angular-momentum transport mecha- 
nisms, including, in particular, radial transport by gravita- 
tional torques and by MRI-induced turbulence. In this paper 
we consider only vertical transport by ccntrifugally driven 
winds in an attempt to model a radially localized disc re- 
gion where this mechanism may dominate. (Note, however, 
that both vertical transport and radial transport notably 
MRI-induced turbulence - could in principle operate at the 
same disc radius; see Salmeron, Konigl & Wardle 2007.) As 
discussed in Paper I, the large-scale, ordered magnetic field 
envisioned in this scenario could be either interstellar field 
advected by the accretion fiow or dynamo-generated field 
produced in cither the star or the disc. In view of the strong 
evidence for strong outfiows from the inner regions of pro- 
tostellar discs, we also neglect alternative modes of angular 
momentum transport that could be mediated by such a field, 
including magnetic braking, 'failed' winds and non-steady 
phenomena. Our treatment has been deliberately simplified 
to facilitate comparison with the analytic results of Paper I; 
in particular, we assume that the matter is everywhere well 
coupled to the field (i.e. A > 1) and that the same conduc- 
tivity regime applies between the mid-plane and the disc's 
surface. In reality, the disc may be weakly coupled between 
the midplane and some finite height and its diffusivity prop- 
erties are expected to change with 2 (e.g. Salmeron & War- 
dle 2005) . Our approximation should be adequate for repre- 
senting the behaviour of the dominant diffusivity regime in 
the well-coupled region of the disc. However, if a significant 
fraction of the local column density is magnetically weakly 
coupled then the vertically averaged properties of the disc 
(such as the infiow speed) could be significantly modified 
(see Li 1996; Wardle 1997). 

In conclusion, the results presented in this paper con- 
firm the validity of the parameter constraints derived in Pa- 
per I for physically viable configurations of Hall diffusivity- 
dominated protostellar discs in which ccntrifugally driven 
winds dominate the local angular momentum transport. 
They also demonstrate that the algebraic expressions de- 
rived on the basis of the hydrostatic approximation correctly 
identify the generic properties of such discs and arc useful for 
clarifying the behaviour of the full numerical solutions. More 
generally, the theoretical framework developed in WK93, 
Paper I and the present work can be used to study discs 
of this type also in other diffusivity regimes and in other as- 
trophysical environments. In particular, it can help interpret 
observations of such systems by relating the properties of the 
outflow to those of the underlying disc (e.g. Konigl 2010). 
It may also be useful for guiding non-ideal-MHD numerical 
simulations of wind-driving discs. 
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APPENDIX A: PARAMETER CONSTRAINTS 
IN THE HALL LIMIT 

The parameter constraints on viable wind-driving disc solu- 
tions were obtained in Paper I by applying the hydrostatic 
approximation to a radially localized disc model under the 
assumption that the charged component of the weakly ion- 
ized disc material consisted of two particle species: positive 
'ions' (subscript 'i') and negative 'electrons' (subscript 'e'), 
each singly charged. Instead of employing two independent 
ratios of the conductivity-tensor components as model pa- 
rameters, as is done in this paper, we used the ion and elec- 
tron Hall parameters, i.e. the ratios of the gyrofrequency and 
the collision frequency of these two species with the neutrals, 
which are given by 

o eB 1 . eB I ^ , . ^ , 

/3e = , A = =qPe, (Al) 

ruic 7ip rUeC 7eP 

where c is the speed of light, m is the particle's mass, e is the 
unit electric charge and the total mass density p is used to 
approximate the neutral mass density. The collisional cou- 
pling coefficient 7j is equal to {av)j/{mj + m), where m 
is the mean mass of the neutral particles and {ov)j is the 
rate coefficient of momentum exchange of species j with the 
neutrals. The product 7ip (7ep) in equation (Al) therefore 
represents the momentum-exchange collision frequency of 
the ions (electrons) with the neutrals. In the above defini- 
tions, B = |B| sgn{Bz}, so the sign of the Hall parameter is 
sensitive to the magnetic field polarity. It is further assumed 
that the ions are 'heavy' and the electrons are 'light', so that 
(?< 1.® 

Another key parameter employed in Paper I (see also 
WK93) is 



T ; 



7iPi 
Ok 



(A2) 



(where pi is the ion mass density) , the ratio of the Keplerian 
rotation time to the neutral-ion momentum exchange time. 
In the ambipolar diffusivity limit (ao ^ ap ^ Ian] or, 
cquivalently, |/?i| S> 1) the Elsasser number A (Section 2.2) 
reduces to T, whereas in the Hall diffusivity limit (up <C 
|o-h| < o-o or, equivalently, |/3i| < K |^e|) A = T|^i|. 

The parameter constraints for the four Hall sub-regimes 
are presented in Table Al, which reproduces Table I.l. The 
physical origin of the imposed constraints is summarized in 
the caption of this table. The key predicted properties of the 
solutions in each sub-regime are listed in Table A2, which 
reproduces Table 1.2. We also reproduce below some of the 
expressions used in the derivation of these results that are 
relevant to the analysis presented in this paper. 

The ratio \dhr/dh4,\o is given by 

\dW\ _ 2To + /? 



1 + . 



® In a real disc containing two charged species, the value of q is 
fixed by the physical properties of the charge carriers (see equa- 
tions 1. 7-1. 9). This, in turn, constrains the values that the con- 
ductivity ratios cth/^x and ctx/ctq can take (see equations 84 
and 85) . We do not incorporate this constraint into the parameter- 
space analysis in Section 4 so as not to unduly complicate the 
discussion; however, we have verified that the 'effective' value of 
q, obtained from the ratio of equation 84 to equation 85, remains 
<g; 1 for all the solutions that we present. 
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Table Al. Parameter constraints for wind-driving disc solutions in the limit where the Hall diffusivity dominates and assuming ee = 0. 
Four distinct cases can be identified, depending on how the values of so = j3eoAo a^id of 2Ao = 2To|/3io| compare with 1. The first 
inequality expresses the requirement that the disc remain sub-Keplerian below the wind zone (z < the second is the wind launching 
condition (the requirement that the magnetic field lines be sufficiently inclined to the vertical for centrifugal acceleration to occur), the 

third ensures that the base of the wind is located above the (magnetically reduced) density scalehcight and the fourth specifies that the 
rate of Joule heating at the midplane should not exceed the rate of release of gravitational potential energy there. 



Case Limits Parameter Constraints - Hall Limit 

So = /3eoAo ^0 = T~o|/3io| (multi-fiuid formulation) 



(i) 
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(iv) 
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Table A2. Key properties of viable disc solutions in the Hall regime. Listed, in order, are the midplane values of \dhr/db^\, the 
magnetically reduced scaleheight in units of the tidal scaleheight {h = h/hT), the similarly normalized vertical location of the base of 
the wind % ™ units of h and the normalized Joule dissipation rate j • e' at the midplane. 



Case 


So = /3eO 


Limits 
fto Ao = To|Ao| 


\dbr/db,j,\o 


Solution Characteristics - Hall Limit 
h zi,/h 
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(equation L115), where /9 = 1/pio- This expression can be 
used to approximate \brh/bri,h\, the ratio of the corresponding 

magnetic field components at the base of the wind. Another 
relationship between the field components at Zb is provided 
by 



(A4) 



involves the ratio of the tidal scaleheight to the disc radius, 
hr/r = Cs/vk- This parameter does not appear explicitly in 

the normalized equations for the disc structure (Section 2.1), 
but it is used in matching the disc solution to a self-similar 
wind solution (via equations 62 and 63). We have verified 
that this constraint is satisfied by the matched disc/ wind 
solutions presented in Section 4.3. 



(equation 1.90). Since \db^/dbr\o is always < 1 in the Hall 
regime, one can approximate, to leading order, brh V^/ao 
and b^h ~ — \/2(l + qj3^)/[{2To + j3)ao\. The assumption 
\db^/dbr\o <^ 1 also leads to the following simplified expres- 
sion for the magnetically reduced scaleheight: 



2ao l + 



2To + P 



(A5) 



(equation 1.119). The height of the base of the wind is, in 
turn, given by 



aot Tg-h(5/2)To/3-|-/3' 



3^/2Ta 



l + g/32 



(A6) 



(see equation L120). Finally, the Joule dissipation term in 
Table A2 is evaluated from 



{j • e')o = 



4Toa§ 



(i + gr) 



1 + 



2To + / 
l + qff- 



(A7) 



(see equation L122). The rate of Joule heating at the disc 
mid-plane should not exceed the rate of release of gravita- 
tional potential energy at that location. This requirement 
is expressed by the rightmost inequality of Table Al, which 
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